Doctor of Philosophy

dissertationPlatelet aggregation, an important part of the development of blood clots, is a complex process involving both mechanical interaction between platelets and blood, and chemical transport on and o the surfaces of those platelets. Radial Basis Function (RBF) interpolation is a meshfree method for the interpolation of multidimensional scattered data, and therefore well-suited for the development of meshfree numerical methods. This dissertation explores the use of RBF interpolation for the simulation of both the chemistry and mechanics of platelet aggregation. We rst develop a parametric RBF representation for closed platelet surfaces represented by scattered nodes in both two and three dimensions. We compare this new RBF model to Fourier models in terms of computational cost and errors in shape representation. We then augment the Immersed Boundary (IB) method, a method for uid-structure interaction, with our RBF geometric model. We apply the resultant method to a simulation of platelet aggregation, and present comparisons against the traditional IB method. We next consider a two-dimensional problem where platelets are suspended in a stationary fluid, with chemical diusion in the fluid and chemical reaction-diusion on platelet surfaces. To tackle the latter, we propose a new method based on RBF-generated nite dierences (RBF-FD) for solving partial dierential equations (PDEs) on surfaces embedded in 2D domains. To robustly tackle the former, we remove a limitation of the Augmented Forcing method (AFM), a method for solving PDEs on domains containing curved objects, using RBF-based symmetric Hermite interpolation. Next, we extend our RBF-FD method to the numerical solution of PDEs on surfaces embedded in 3D domains, proposing a new method of stabilizing RBF-FD discretizations on surfaces. We perform convergence studies and present applications motivated by biology. We conclude with a summary of the thesis research and present an overview of future research directions, including spectrally-accurate projection methods, an extension of the Regularized Stokeslet method, RBF-FD for variable-coecient diusion, and boundary conditions for RBF-FD

A High-Order Radial Basis Function (RBF) Leray Projection Method for the Solution of the Incompressible Unsteady Stokes Equations

A new projection method based on radial basis functions (RBFs) is presented for discretizing the incompressible unsteady Stokes equations in irregular geometries. The novelty of the method comes from the application of a new technique for computing the Leray-Helmholtz projection of a vector field using generalized interpolation with divergence-free and curl-free RBFs. Unlike traditional projection methods, this new method enables matching both tangential and normal components of divergence-free vector fields on the domain boundary. This allows incompressibility of the velocity field to be enforced without any time-splitting or pressure boundary conditions. Spatial derivatives are approximated using collocation with global RBFs so that the method only requires samples of the field at (possibly scattered) nodes over the domain. Numerical results are presented demonstrating high-order convergence in both space (between 5th and 6th order) and time (up to 4th order) for some model problems in two dimensional irregular geometries.Comment: 34 pages, 8 figure

Effects of oxygenate and aromatic content on engine-out aldehyde emissions from pure, binary, and ternary mixtures of ethanol, toluene, and iso-octane

Sustainable fuel components, such as ethanol, can be blended into gasoline to help reduce fossil fuel consumption. Aldehydes are both observed emissions and major intermediates in the oxidation of gasoline/ethanol fuel mixtures and are solely attributed to the fuel’s alcohol content. This study aims to evaluate the direct impact of toluene, iso-octane, and ethanol on engine-out formaldehyde and acetalde hyde emissions. A single-cylinder direct injection spark ignition engine was run at low speed and load conditions with varying equivalence ratios. The emissions are measured using a FTIR. The aldehyde emissions from pure ethanol are ten- and five- times greater than pure toluene and iso-octane, respectively. Greater formaldehyde than acetaldehyde is found for iso-octane and toluene and vice versa for ethanol. The addition of 25 %vol of toluene to ethanol halves the average aldehyde emissions due to toluene suppressing effects. In the ternary mixtures with fixed ethanol content, the higher toluene concentration mixture led to lower average aldehyde emissions

Mesh-Free Semi-Lagrangian Methods for Transport on a Sphere Using Radial Basis Functions

We present three new semi-Lagrangian methods based on radial basis function (RBF) interpolation for numerically simulating transport on a sphere. The methods are mesh-free and are formulated entirely in Cartesian coordinates, thus avoiding any irregular clustering of nodes at artificial boundaries on the sphere and naturally bypassing any apparent artificial singularities associated with surface-based coordinate systems. For problems involving tracer transport in a given velocity field, the semi-Lagrangian framework allows these new methods to avoid the use of any stabilization terms (such as hyperviscosity) during time-integration, thus reducing the number of parameters that have to be tuned. The three new methods are based on interpolation using 1) global RBFs, 2) local RBF stencils, and 3) RBF partition of unity. For the latter two of these methods, we find that it is crucial to include some low degree spherical harmonics in the interpolants. Standard test cases consisting of solid body rotation and deformational flow are used to compare and contrast the methods in terms of their accuracy, efficiency, conservation properties, and dissipation/dispersion errors. For global RBFs, spectral spatial convergence is observed for smooth solutions on quasi-uniform nodes, while high-order accuracy is observed for the local RBF stencil and partition of unity approaches